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ABSTRACT 

We measure the UV-optical color dependence of galaxy clustering in the local universe. Using 
the clean separation of the red and blue sequences made possible by the NUV — r color- 
magnitude diagram, we segregate the galaxies into red, blue and intermediate "green" classes. 
We explore the clustering as a function of this segregation by removing the dependence on 
luminosity and by excluding edge-on galaxies as a means of a non-model dependent veto of 
highly extincted galaxies. We find that ^(rp, tt) for both red and green galaxies shows strong 
redshift space distortion on small scales - the "finger-of-God" effect, with green galaxies 
having a lower amplitude than is seen for the red sequence, and the blue sequence showing 
almost no distortion. On large scales, ^(rp, tt) for all three samples show the effect of large- 
scale streaming from coherent infall. On scales \h~^ Mpc < < lQh~^ Mpc, the projected 
auto-correlation function Wp{rp) for red and green galaxies fits a power-law with slope 7 ^ 
1.93 and amplitude tq ^ 7.5 and 5.3, compared with 7 ^ 1.75 and tq ^ 3.9h^^ Mpc for 
blue sequence galaxies. Compared to the clustering of a fiducial L* galaxy, the red, green, 
and blue have a relative bias of 1.5, 1.1, and 0.9 respectively. The Wp{rp) for blue galaxies 
display an increase in convexity at ^ Ih^^ Mpc, with an excess of large scale clustering. Our 
results suggest that the majority of blue galaxies are likely central galaxies in less massive 
halos, while red and green galaxies have larger satellite fractions, and preferentially reside in 
viriaUzed structures. If blue sequence galaxies migrate to the red sequence via processes like 
mergers or quenching that take them through the green valley, such a transformation may be 
accompanied by a change in environment in addition to any change in luminosity and color 

Key words: methods: statistical - galaxies: elUptical and lenticular - galaxies: evolution - 
galaxies: clusters: general 
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1 INTRODUCTION 



With the advent of the Sloan Digital Sky Survey fSDSS lYorket alJ 
l200Ch and its value added galaxy catalogs, it has been possible to 
study the subject of galaxy bimodality and its relationship to fun- 
dam ental properties, such as stellar mass and star-formation history 
(e.g.,lKauffmann et al.l2004ISchiminovich et al.l2007l : [Salim et alJ 
l2007h . The broad division of galaxies into star forming disks and 
quiescent early type galaxies is the fundamental principle of Hub- 
ble's tuning fork system of classification and is well established. 
In a plot of optical g — r color vs Mr, red galaxies define a clear 
sequence, while the locus of blue galaxies is broadened into the so- 
called "blue cloud". The red seq uence has been shown to maintain 
its integrity with look-back time jBower et aljl 992) and h as grown 
in mas s since redshift ^ 1. St udies by Bell et al. ( 2004h . lBlantoj 
1 I2OO6I) , iFaber et all lioOTh . and lBrown et al., C2008.) argue that the 
stellar mass contained within the red population has increased by 
roughly a factor of two in half the Hubble time. 

A significant breakthrough in expressing this blue/red di- 
chotomy occurred when photometry from the Galaxy Evolu- 
tion Explorer (GALEX), notably the near-ultra violet ( NUV) 
band, was matched with SPSS photo metry (Martin etal. II2OO7I : 
IWvder et alj |2007| : ISchiminovich et al.l l2007l : ISalim et aTmOOTh . 
When the diagram is plotted using NUV — r as the color, two 
clear sequences emerge: the familiar red sequence and a new blue 
sequence in place of the "blue cloud" of optical studies. Between 
the blue and red sequences there are galaxies present in a so- 
called "green valley". Many of thes e are spectroscopically classi- 
fied Type II active gala ctic nuclei dRich et al.ll2005l : iMartin et al.l 
I2007l;fsalimetalj2007l). 

Ipaber et al.' ('2007t), IMartin et al.l j2007h . and 

ISchimino vich et al. ( 2003) propose several paths by which 
galaxies might transition from the blue to the red sequence. The 
presence of AGN in the green valley suggests that AGN activity 
is associated with a quenching of star formation jSilk & Reed 
ll998l : lHopkins et alj|200d 2007). Other paths from the blue to the 
red sequence might, hypothetically, inv olve gas-rich mergers of 
blue galaxies jToomre & ToomrdI 1 1972), or v irial shock heating 
of cold gas streams /Pe kel & Bimboinj||200d) . The red sequence 
might consolidate in luminosity via dissipationless mergers of red 
galaxies, or low luminosity blue galaxies might acquire bulges 
through mergers with starbursts, retaining sufficient mass to land 
the evolved galaxy on the red sequence. However, the green valley 
might also be populated by casual visitors — red galaxies that 
acquire gas and form stars or feed a central AGN. Following this 
brief burst of star formation, these galaxies might ultimately return 
to the red sequence from which they started. 

In lSaiim et al.l l2007h . a plot of mass against specific star for- 
mation rate reveals a clear division between lower mass, star form- 
ing, blue sequence galaxies, and more massive AGN, which are not 
detected in large numbers until stellar mass M ~ 3 X 10" Mq. 
The process responsible for populating the green valley and for po- 
tentially contributing to evolution from blue to red must bear some 
relationship to environment and to the dark matter halos in which 
the galaxies reside. In this study, we investigate the clustering en- 
vironment of the blue and red sequences, and for the green valley. 

The current paradigm of structure formation assumes that 
galaxies are assembled in dark-matter halos. The dependence of 
the clustering on galaxy properties may provide clues to the bary- 
onic processes that are important to galaxy formation and evo- 
lution. The dependence of galaxy clustering on galaxy type has 
been known since the earliest studies of extragalactic astronomy 



( lHubble|[l933 ; IZwickv et"^ 19681). I n the modern era of large- 
scale galaxy surveys, |Pavis&Gelle3 ( ll976l) showed that the an- 
gular auto-correlation of ellipticals has a steeper power-law slope 
that those of spirals. Recent redshift surveys using the Two- 
Pegree Field Galaxy Survey (2dF) and SPSS confirms these ear- 
lier results and t he apparent bimodal nature of galaxy clustering 
kladgwick et al ] I2OO3I: iBudavari et al.1 l2003l : IZehavi et al.1 l200a - 
Iti et alJI2006ij :iWang et al.' 2007h . 

Studies using SPSS have further revealed tha t galaxy color 
is the p roperty most predictive of local environment. iBlanton et al.l 
( l2005t j) found that at fixed luminosity and color, density does not 
correlate with surface brightness nor the Sersic index, and argue 
that morphological properties of galaxies are less closely related 
to environment than their star- formati on history, and are traced by 
broadband optical colors. (See Park et al. (2007) for an alternative 
analysis and point of view.) iLi et alj 12006 a) found that the de- 
pendence of clustering on optical g — r color and D4000 is much 
stronger than structural parameters like concentration and surface 
brightness, and extend to 5/i~^ Mpc, beyond what is expected from 
the localized halo paradigm of structure formation. They concluded 
that at fixed stellar mass, the clustering properties of the surround- 
ing dark matter haloes are correlated with the color of the selected 
galaxies. They further argued that different physical processes may 
be required to explain environmental trends in star formation, dis- 
tinct from those established by galaxy structure. 

In this paper, we will consider the color dependence of the 
two-point physical correlation function of galaxies using samples 
constructed from GALEX, augmented with redshift and optical 
data from SPSS. In particular, we use the natural separation from 
the NUV — r color to assign galaxies into three subsamples of 
red, green and b lue galaxies. Our study complements recent work 
by H einis et"aD (|2009) who investigate the physical clustering of 
galaxies as a function of star- fo rmation history in the lo cal universe, 
as well as earlier studies b y iMilliard et al ] ( l2007l) . I Heinis et al.l 



(2007) and iBasu-Zvch etai] ( l2008h who investigate the angular- 
correlation function of rest frame C/V-selected galaxies and their 
evolution. We measure the auto-correlation function each of the 
different subsamples of galaxies, as well as the cross-correlation 
function between the subsamples. In section|2l we describe in detail 
the data used in this analysis. In section |3| we describe the method 
for estimating correlation functions. We present our results in sec- 
tion|4l and discuss their implications for the nature of green valley 
galaxies and the formation of red sequence galaxies. We summarize 
our findings in section |6| 



2 DATA 

2.1 GALEX and SDSS Data 

The ultraviolet imaging portion of the data-set is from the 
Galaxy Evolution Explorer Satellite (GALEX) that was launched 
in 2003 April (Martin etal. 2005; Morrissev et al. 2005, 200^). 
GALEX obtains wide field imaging in both the far-UV (FUV; 
centered at 1540A) and the near-UV (NUV, 2300A) over a 1.2° 
diameter field of view, with 5" images. Here we use data from 
the Medium Imaging Survey (MIS); these images are ft; 1500 
sec in duration reaching NUV ft; 23 mag, covering an orbital 
shadow crossing. The MIS pointing that defines our sample tar- 
gets the North Galactic Cap, which overlaps the SPSS spectro- 
scopic footprint; this part of the program was designed from stud- 
ies cross-matching SPSS and GALEX data. The data-set used for 
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our current analysis is from the Galaxy Release 3 (GR3) which is 
available from the Multi-mission Archive at STScI (MAS T). The 
GALEX pipeline uses SExtrator teertin & Amoutj 1 19961) to de- 
tect sources and measure fluxes. We use the "MAG_AUTO" output 
from SExtra tor as our default flux measurement; it is essentially a 
lKronl ( ll98(]h magnitude with an elliptical aperture. 

Because our analysis requires each galaxy to have a spectro- 
scopic redshift, we start our cross-matching with a galaxy from the 
SDSS Main spectroscopic survey. For each SDSS galaxy, we search 
for the closest GALEX detection within 4" radius from the loca- 
tion of the SDSS spectroscopic fiber Only GALEXsources within 
0.55° of the tile center field-of-view (FOV) are retained, since as- 
trometry degrades toward the periphery of the FOV (a problem 
which will be resolved in later releases) and the incidence of ar- 
tifacts increases as well. After the matching of GALEX and SDSS 
sources, we further trim the sample to create a statistically complete 
data-set following the procedure of Wyder et al. (2007). We will re- 
fer the reader to their table 1 for the full details. Here, we list a few 
essential parameters and the minor modifications we employed: 
14.0 < r < 17.6, 0.03 < z < 0.25, a,- < 0.2, Zconf > 0.67, 
GALEXexposure time t > 750s and 16.0 < NUV < 23.00 
All magnitudes are AB magnitudes and corrected for Galactic fore- 
ground extinction. 



2.2 SDSS Large-Scale Structure Sample 

Considerable effort has been invested by the SDSS team to prepare 
the redshift data for large-scale structure studies. This secondary 
data-set known as the New York University Va l ue Add ed Catalog 
(NYU-VAGC) is documented in lBlanton et alj j2005al) and avail- 
able from the NYU website. Proprietary versions of this catalogue 
have been used by many groups within the SDSS collaboration 
for various investigations of clustering and luminosity function of 
galaxies. We match the GALEX-SDSS catalog constructed above 
with the SDSS large-scale structure (LSS) DR5 sample. The ver- 
sion used for our analysis includes all of the detailed radial and 
angular selection functions fo r the various statis tical subsamples 
used in previous analyses (e.g. IZehavi et al.ll2005l) . T his sample is 
ideal b ecause we can compare our result with that of IZehavi et al.l 
( I2OO5I) which was based solely on selection from optical criteria. 



2.3 GALEX-SDSS Overlapping Footprint 

In order to statistically define our combined GALEX-SDSS sur- 
vey, we need to have the understanding of the angular sampling 
function of the two surveys, which varies across different regions of 
the sky. GALEX' s Medium Imaging Survey (MIS) survey consists 
of overlapping circular tiles (radius 0.55°), while the SDSS spec- 
troscopic survey is a combination of circular spectroscopic plates, 
but with fiber placement based on a rectangular imaging survey that 
runs along great circles. To combine the footprints of both surveys, 
we use the GESTALT footprint server. GESTALT uses a hier- 
archical pixelization system, enabling one to encode observations 
of arbitrary geometry while tagging information about complete- 
ness and the masking of artifacts. We first encode the GALEXMIS 
survey using GESTALT. We then obtain the detailed observational 



Our faint-end limit of 23.0 is about 0.5 mag fainter than those employed 
bv lWvderet ai] j2007). 

^ http://nvogre.phyast.pitt.edu: 8080/gestalt_tutorial/ 



footprint of the SDSS Large-scale structure sample from the NYU- 
VAGC website. The SDSS footprint is exp ressed as a set of dis- 
joint polygons using the software MANGLE ( [Hamilton & TegmarkI 
2002) which takes into account the complex angular mask and ge- 
ometry of the SDSS survey. We convert these polygons into the 
pixelization scheme of GESTALT and consider the intersection 
of the two surveys. There are approximately 490 square deg in 
the GALEX-SDSS overlapping footprints after masking for holes, 
bright stars and satellite trails, and excluding defects. 

In order to measure the correlation function, a random sample 
needs to be constructed to normalize the galaxy pair counts. Angu- 
lar sampling completeness as a function of position in the sky en- 
coded in the footprint server is used to generate random samples of 
density roug hly 50 times the galaxy density. We adopt the method 
proposed by iLi et al.l ( l2006ah where we assign each galaxy in our 
sample to a random position on the sky but keep all other attributes 
the same (e.g. redshift, color, magnitude). This random sample by 
construction has the same redshift distribution as the original sam- 
ple, and thus does not smooth out the redshift st ructure like thos e 
generated via the luminosity function. As noted in iLietaljjiooeah . 
this approach works well in surveys with a wide-angular sky cov- 
erage (e.g. much larger than the typical large-scale structure), and 
with small variation in survey depth. We note here that our random 
sample would inherit the redshift correlation function of the color- 
magnitude distribution of the parent galaxy sample, only spatial 
distribution has been randomized, hence any excess in clustering 
must be due to positional differences. 

In the SDSS spectroscopic survey, no two galaxies with sep- 
aration 6 less than 55" can both be assigned spectroscopic fibers 
for observation on any given observing plate. Hence, a large frac- 
tion of galaxy pairs with 8 < 55" are missing. We correct for this 
"fiber-collision" problem by using the observed arigular correlation 
function, a method first suggested [Li et al. ( 2006a p, and described 
in detail in the companion paper by Heinis et al.l ( l2009h . In brief, 
the observed projected two-point angular correlation function for 
both the photometric sample lUph(^) and the spectroscopic sample 
Wsp (8) are measured and used to construct the pair weighting ratio: 

FiO) = ^^^Hllgl±i (1) 
Wsp(8) + 1 

as a function of separation. Empirically, iHeinis et al.l ('2009^ finds 
F{e) ~ 3 for e < 55" and zero otherwise for the all GALEX- 
SDSS cross-matched galaxies. We apply this correction to each of 
the red, green and blue subsamples equally. 



3 METHODOLOGY 

3.1 The Two-point Correlation Function 

In brief, the two-point auto-correlation function ^(r) measures the 
excess probability of finding a galaxy pair with separation r from a 
random galaxy distribution, 

dP = n[l + ^{r)]dV (2) 

where n is the mean number density of galaxy sample. For the 
last forty years, the correlation function has served as the pri- 
mary method for cosmologists to quantify the clustering proper- 
ties of galaxie s from large scale surveys (Totsuji&Kihara 1961; 
|Peebleslll980l) . If the underlying density distribution is Gaussian, 
then the correlation function fully describes all statistical properties 
of a given distribution. Since early cosmological models are often 
based on primordial Gaussian dark matter density fields, the use 
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of the correlation function leads to a straightforward comparison 
between empirical studies on the statistical distribution of galax- 
ies with such models. Recently, with the wide spread adoption of 
a standard A-dominated cosmology I Ostriker & Steinhardt 1995; 
iRiess et alJ[T993 : IPerlmutter et alj|l99'9l ; ISpergel et allboOTl) . the 
correlation function is used instead to probe the growth of structure 
in the universe and the range of formation scenarios for varying 
kinds of galaxies. The correlation length - the amplitude from a 
power-law correlation function - provides information on the mass 
of dark-matter halos in which the various galaxies reside, link- 
ing observation wi t h theoretical descri ption of structure formation 
( iBower et al.ll2006l : ICroton et alj2006l) 

For a galaxy survey with a well-defined angular selection 
function, ^ can be esti mated using an optimal estimator like the 
iLandv & Szalavl ( Il993h estimator: 



£,LS 



DD - 2DR + RR 
RR 



(3) 



where DD, DR and RR are normalized counts of galaxy pairs 
in the data-data, data-ra ndom and random-random catalogs|j The 
iLandv & Szalavl l l 19931) estimator is preferred because it is rela- 
tively in sensitive to the size of the random catalog and to edge cor- 
rections jKerscher et alj2OO0l) . 

Because we observe galaxies in redshift and projected space 
and not in physical space, in practice, the correlation function is 
first measured on a two-dimensional grid of separations: n along 
the line of sight (redshift space) and rp for angular separation on 
the sky. In addition to providing information about the underly- 
ing mass distribution through the amplitude of the clustering sig- 
nals, the two-dimensional correlation function ^(7r,rp), contains 
additi onal information about the dynamics of the galaxies jPeeblesI 
119801) . At small projected separations rp, random motions within 
virialized over-density (e.g. clusters of galaxies) causes an elonga- 
tion along the line-of-sight (vr direction) known as the "finger-of- 
God" effect. On large scales, coherent streaming of galaxies into 
potential wells c auses an apparent compr e ssion of structure along 
the li ne-of-sight l lSargent&Tumedll977l : lKaisedll987l : iHamiltoiJ 
Il99 j ). Various studies have use d ^(tt, rp) to extract cosmological 
and dynamical information (e.g., lTinkeilf2007l) . 



3.1.1 Real-Space Correlation 

Because redshift space distortion only affects the line-of-sight com- 
ponent of ^(tt, rp), we can recover the true space correlation func- 
tion ^(r) by following the standard procedure of computing the 
projected correlation function: 



poo 

Wp{rp) = 2 / d-n(_{rp,-n). 
Jo 



(4) 



In practice, we integrate along the line-of-sight direction out to 
TT = 30h~^ Mpc. This is large enough to include almost all cor- 
related pairs but also stable enough to suppress noise from distant 
uncorrelated pairs. The projected correlation function can in turn 
be related to the real space correlation function ^(r). 



2 

Wp{rp) = 2 / rdr^{r){r — rp) 
Jo 
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(5) 



If the real-space correlation function follows a power-law ^(r) = 
(r/ro)^, we can infer its parameters: the correlation length ro and 
the power-law slope 7 from the best-fit power-law to Wp{rp) using 
the following deprojection: 



■Wp{rp) =rp(-^ 



(6) 



3.1.2 Cross-correlation 

Related to the auto-correlation function is the cross-correlation 
function between two classes of galaxies. The cross-correlation 
function ^1,2(7'), measures the clustering of one type of galaxy 
around another. ^1,2 (r) is essentially the probability of finding a 
galaxies of type 1 around a galaxy of type 2 as a function of sep- 
aration r. F or our analysis, we use t he cross-correlation version of 
the classical lOavis & PeeblesI ( Il983h estimator: 

? - ^'^^ 1 m 

^'''-Djh'^ 

where D1D2 are the normalized counts of cross-pairs; while D1R2 
are cross-pairs of type 1 with random galaxies having the same 
redshift distributions as galaxies of type 2. 

Recently, the cross-correlation function h as been used exten- 
sively in clustering studies of galaxy properties jZ ehavi et al. 200^; 
|Wang et alj|2007l : ]coil et"ai]|2008l ; IPadmanabhan et al...2008 ; .Che3 
2009[) since the auto-correlation function alone tells us little about 
how galaxies of various types relate to one another. Two popula- 
tions of galaxies may both have comparable correlation strength, 
and yet be physically unrelated if they are spatially segregated. 
Specifically, the cross-correlation function may be used in conjunc- 
tion with the auto-correlation functions of the galaxies to glean 
information as to how they mix statistically. For example, if the 
two populations are well mixed, i.e. they are consistent with being 
drawn stochastically from their respective auto-correlation function 
equally, their projected cross-correlation function Wp''^{rp) would 
follow that of the geometric mean of their auto-correlation func- 
tions 



1,2/ 

Wp' [rp 



Wp^{rp)wl'^{rp) 



(8) 



On the other hand, if the populations were segregation between the 
populations and do not distribute evenly in all space, e.g. they sit 
on different halos, the amplitude of Wp''^{rp) would be lower than 
that of Wp''^{rp) known in the literature as stochastic anti-bias (e.g. 
iBlantoa.200Q ; .Swanson et al..,2008i) 

3.1.3 Bootstrap Errors 

We estimate the er rors for the correlation measurements using a 
modified bootstrap tefronll98lh method for spatial statistics known 
as marked-point bootstrap ( lLolj2 008). We first estimate the contri- 
bution of each galaxy to the correlation function - the marks: 

n 

^^('■) = 5Z (t}{\xi~Xj\e{r-dr,r + dr)} (9) 
such that 



^ The use of congelation function to probe the growth of stnicture predate 
the advent of A-CDM, especially in the study of faint blue galaxies (e.g. 



lEfstathiouet all 19 91). 
* Normalized counts of galaxy pairs are weighted by the selection function 
of the galaxies involved. 



(10) 



where n is the number of galaxies and (j> is the chosen correlation 
estimator (e.g. ^ls)- Hence, the mark associated with galaxy Xi 
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is the number of excess pairs at a distant r from Xi. These are 
(roughly) independent, and identically distributed, and can be used 
to replicate B bootstrap samples by the sampling with replacement. 
An estimate of the correlation function ^* can be computed from 
each of the j-th bootstrap samples. We use B — 999 bootstrap 
samples for our analysis. The distribution of £,*{ri) can be use to 
estimate the error bars of the correlation function at each separation 
ri. For an equivalent of one sigma errors, we ranked ^* (r,) and take 
the 159"'' and 840"" for the lower and upper error bounds. Because 
the ^ errors for the are correlated, when we estimate parameters 
for the two parameter power-law model (equation we refitted 
each bootstrap sample to obtain a two dimension distribution of the 
parameters. 

We chose this method over the conventional jackknife meth- 
ods because of the fragmented nature of the GALEX-SDSS foot- 
print. Our bootstrap errors are consistent with the jackknife er- 
rors estimated by excluding one GALEX tile at a time. We note 
here that jackknife errors are known to overestimate the vari- 
ance o n small-scale, and often bias the overall correlation esti- 
mates jNorberg et alj2009l) . Since bootstrap is a form of internally 
estimated errors, and internally estimated errors are known not 
to rep roduce externally estimated errors accurately (Norberg et al. 



to rep : 
\2Q0% 



our clustering results should be treated with caution. 



4 DEPENDENCE OF CLUSTERING ON (7y-OPTICAL 
COLOR 

Figure [T] shows the NUV — r color-magnitude diagram (CMD) in 
the local universe. All absolut e magnitudes are fc-correc ted using 
the K_CORRECT program by Isianton & RoweisI ( l2007h and are 
quoted in units with Hubble constant Hq = lOO/i km Mpc~^ 
, i.e. Mr = Mr — 5 log(/i = 1). Compared with the optical color- 
magnitude diagram (e.g. g — r vs Mr), the NUV ~ r color axis 
has a higher dispersive power. The scatter plot shows two clear 
sequences: red and blue, and an intermediate "valley" of "green" 
population. Conventional optical diagrams only display a single 
red sequence with an extended blue cloud. Hence, diagnostic di- 
agrams like Figure[T]separate galaxies into three natural groupings: 
(1) a red sequence of bulge-dominated galaxies with old stellar 
systems, (2) a blue sequence of star-forming systems consisting 
of mainly late-type galaxies, and (3) galaxies in a " green valley" 
that in principle might exhibit t ransitional properties jMartin et al.l 
l2007l : [Schiminovich et alj2007l) . Detailed luminosity functions and 
physical properties analysis of these galaxies derived from this 
C/y-optical CMD has been reported e l sewhere jWvder et al .120071 : 
iMartin et alj |2007| : ISalim et alj l2007l : ISchiminovich et al.l 120071) . 
Here, we investigate the clustering properties of these three sub- 
populations of galaxies, separated by the two tilted horizontal lines, 
using the two-point auto-correlation function, as well as their rela- 
tion to each other using the cross-correlation function. Our results 
are not sensitive to 0.25 mag shift in the oblique equation used to 
classify the galaxies. 



4.1 Dust 

Dust content within each galaxy can modify their broadband col- 
ors, hence the resultant col or-magnitude diagram, as galaxies move 
from one group to another ^Martin et alJI 20071 : ISchiminovich et al] 
l2007h . This is a concern, as it might affect our estimate of the cor- 
relation function by either increasing the scatter between groups, 
or bias the result in a systematic but unknown manner. Indeed, a 



fraction of the green valley galaxies are dusty star forming galax- 
ies whose intrinsic color would have placed them on the blue se- 
quence in the absence of dust l iWvder et al. 1 12007'). However, the 
available procedures for dust correction are highly uncertain. They 
give inconsistent results d epending on whether one uses a primar - 



ily photometric approach dSalim et al.ll2007l: I Johnson et al 



or one based on spectroscopic indices ( iKauffmarm et al 



200e), 



20041) 



The former, for example, has the side effect of reduc ing the red 
sequence number counts substantially jSchiminovich et al. 20 071; 
iHeinis et al]|2009l) . while the latter approach requires spectroscopy 
of modest S/N, which is not available for part of our sample. A 
comprehensive analysis of the clustering as a function of star- 
formation history wi th a dust-corre cted color-magnit ude diagram 
using the method of I Johnson et alj I2OO6I is done in iHeinis et al.l 
( l2009l) . 

Here, we adopt a geometric approach. Because dust lanes 
in galaxies usually appear when viewed edge-on, we can reduce 
the effect of dust on galaxy color by restricting our analysis to 
galaxies that are primarily face-on. The middle panel of Figure [T] 
shows the "dust-corrected" color-magnitude diagram obtained by 
including only galaxies with isophotal (minor-over-major) axis ra- 
tio b/a > 0.6. To the extent that each of the subpopulations have 
the same intrinsic distribution of b/a, removing edge-on galax- 
ies will give a correct mixture of galaxies that mimic the intrinsic 
color-magnitude distribution. Comparing the dust-corrected CMD 
with the CMD on the left panel (hereafter the full distribution), 
the density of green galaxies around Mr ~ —21 is reduced sub- 
stantially, with the contours showing a more pronounced "valley" 
separating the sequence of galaxies. 

The righthand panel of Figure [T] shows the median ratio of 
edge-on (b/a < 0.4) to face-on (b/a > 0.8) as a function of color 
and magnitude, with the same density contours from the middle 
panel. The region with the highest edge-on/face-on ratio falls in 
the green valley and towards the faint end of the luminosity dis- 
tribution jChoi et alj|2007l : iMartin et alj|2007l : [ Schiminovich et al.l 
I2OO7I) . For the galaxies with the magnitude range with —23.5 < 
Mr < —19.0 used to construct our flux-limited sample, the me- 
dian b/a for red and blue sequence galaxies are 0.7 3 and 0.72 
respectively, while for green valley galaxies it is 0.66. IChoi et all 
12007) find that edge-on galaxies are also statistically fainter due 
to the internal extinction, and it affects morphologically late-type 
galaxies (classified by eye) more than early-type galaxies. Hence 
by limiting our analysis to face-on galaxies, we would reduce bi- 
ases associated with differential dimming in addition to the biases 
from color shifts. Following Choi et al. ( 2007), we restrict our anal- 
ysis to galaxies with b/a > 0.6 as a proxy for a galaxy distribution 
derived from a dust-corrected CMD. 



4.2 Luminosity bias 

The clustering of ga l axies is lu minosity depend ent. iNorberg et al.l 
l l2002l) . IZehavi et alj bOOSt) and lLi et alj ( l2006j) show that the am- 
plitude of the projected correlation function Wp increases mono- 
tonically as a function of luminosity (or stellar mass) on scales 
ranging from 0.2h^^ Mpc to lOh^^ Mpc. In order to study the 
color dependence on clustering, or to compare the clustering in- 
trinsic to the membership of discrete subpopulations, it is vital to 
remove this known luminosity dependence. The top two panels of 
Figure |2] shows the luminosity distribution of the red, green and 
blue population of galaxies with Mr < — 19 from our flux-limited 
sample. The series of dashed histograms (in red, green and blue) 
show the original luminosity distributions. Blue galaxies are more 
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Figure 1. NUV — r vs Mr color-magnitude diagram (CMD) of galaxies in the local Universe. The left panel shows the full distribution. The middle panel 
shows the distribution drawn from face-on galaxies with b/a > 0.6, used to correct for the effect of dust on the CMD. The light panel shows the distribution 
of edge-on/face-on ratio overlay onto the CMD contours. The green valley is dominated by region of the CMD with high relative fraction of edge-on galaxies. 
By limiting ourselve to galaxies with Mr < —20.3 (dashed vertical lines) for the volume-limited sample, we also eliminate region of the CMD with the 
highest fraction of edge-on galaxies. 



numerous and less luminous, on average, compare to red and green 
galaxies, reflecting a steeper blue luminosity function at the faint- 
end (Wyder et al. 2007). To remove this luminosity dependence, we 
resample the luminosity histograms to match the number counts 
of the smallest of the three distributions at each magnitude bin. 
The result of this resampling is shown by the solid histograms. 
The left panel shows histograms from the full distribution, while 
the right shows those from a dust-corrected distribution. There are 
4177 (full) and 3148 (dust-corrected) galaxies in each of the resam- 
pled luminosity distribution of red, green and blue galaxies, with a 
common median luminosity of Mr ~ —21.0 , about half a magni - 
tude brighter than Al*, the typical luminosity dSlanton et alj2003h . 
Note that if the fundamental attribute that drive clustering is stellar 
mass, removing the luminosity dependence like what we have done 
here would still leave residual clustering due to the difference in 
mass-to-light ratio of the respective color selected sample. 

4.3 Volume-limited sample 

While the resampled red, green and blue galaxies are matched in 
luminosity, they are not matched in volume and have different 
redshift distributions. This makes the interpretation of the phys- 
ical correlation function problematic. To this end, we select the 
largest possible volume within our catalog, using a redshift cut of 
0.03 < 2 < 0.12 and an additional luminosity cut at Mr < —20.3 
to construct a volume-limited sample[f| The luminosity distribu- 
tions of this sample are shown on the bottom row of Figure |2] As 
before, the left panel is for the full CMD, while the right panel are 
for the dust-corrected (restricted) CMD. Similar to the flux-limited 
case, the (original) blue dashed histograms are substantially differ- 
ent the red and the green, and are weighted more heavily towards 
the faint-end. We resample the histograms to match in luminosity. 
Note that in the case, the procedure outlined in section |42] essen- 
tially amounts to looking for an optimal function W{L) such that 
W {L)${L) gives the minimum luminosity function of the three 

^ This is not strictly volume-limited since NUV is not complete at these 
re dshifts. An addition al weighting is applied using the luminosity function 
of lWvderet ai]i2007l) . 



subsamples. The new red, green and blue volume-limited subsam- 
ples each consists of 1971 (full) and 1226 (dust-corrected) galaxies 
with a common median luminosity of Air ~ 20.9, almost identical 
to the flux-limited case. 



4.4 Red, Green, Blue Auto-correlation 

Figure [3] and |4] shows the two-dimensional auto-coiTelation func- 
tion ^(7r,r'p) as a function of line-of-sight (tt) and projected (r-p) 
separation for the flux-limited and volume-limited samples. The 
panels on the first row in each figure are from the full distribu- 
tion, while the second row are from the dust-corrected distribution. 
The panels are for red sequence (left), green valley (middle) and 
blue sequence (right) galaxies. We have binned n and r-p linearly 
at lh~^ Mpc. For all panels, the grey scale has the same range, 
and the contours are boxcar smoothed at 2h^^ Mpc. The ^(vr, Vp) 
contours indicate the constant probability of finding pairs at a given 
TT and Tp. The heavy solid line marks ^{n, Vp) = 1, and contours 
are spaced with increment (inner) and decrement (outer) by a fac- 
tor of 2. The effect of redshift-space distortion is clearly seen in 
all panels, manifested by their departure from isotropy (concentric 
circles). At small Vp, the contours are elongated along the line-of- 
sight (tt) due to virial motions of galaxies in clusters. At large Vp, 
the contours are compressed in the vr direction due to the coher- 
ent large-scale stream ing as galaxy infall into potential well. We 
recover the results of IZehavi et al.l ( |2005|) : red sequence galaxies 
show the strongest finger-of-God effect and the larger correlation 
amplitude; and all three subpopulations (of all samples and distri- 
bution) show clear signs of large-scale compression. 

We see a clear finger-of-God effect for green valley galaxies, 
but not for blue sequence galaxies. In Figure[3] if one compares the 
first contour inwards from ^ = 1 (heavy line) for red galaxies to 
the ^ = 1 contour of green galaxies, green and red galaxies have 
identical kinematics, differing only by a scaling in the amplitude. 
The blue sequence appears to have a dynamical structure domi- 
nated by large-scale streaming distinct from both r ed and green 
galaxi es, similar to the optical blue cloud sample of IZehavi et al.l 
l|2005t). We will consider the implications of this in the discussion 
below. The same conclusion can be drawn from the volume-limited 
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Figure 2. The luminosity distribution of our flux-limited (top row) and volume-limited (bottom row) samples. The two panels on the left are for the full 
distribution, the right is for dust-corrected CMD, using only galaxies with b/a > 0.6. The dashed histograms show the original luminosity distribution for 
each of the subsamples of red, green and blue galaxies. We resampled the luminosity to match the number counts of the minimal of the three distribution at 
each magnitude bin, shown as solid histograms. Even for the volume-limited samples, the variations in the original luminosity distribution among the three 
subsamples are substantial. Blue galaxies are heavily weighted towards the less luminous end compared to red or green galaxies. 



sample shown on in Figure |4] The dust-corrected samples (bottom 
row) is substantially noisier We also note that the contours of the 
blue sequence of the volume-limited full sample (top right panel) 
are roughly circular for ~ 4 < iv,rp < lOh^^ Mpc where the 
finger-of-God elongation is balanced by compression due in infall. 

Figure |5] and |6] shows the projected correlation function 
Wp{rp) for the flux-limited and volume-limited samples. The col- 
ored (red, green, blue) points are the respective measured corre- 
lation function. The panel on the left (right) is for the full (dust- 
corrected) distribution. As expected, in all four panels, the red sam- 
ple has the largest amplitude, the blue the lowest, and green in be- 
tween. On large scales (r^ > lh~^ Mpc), the ratio of projected 
auto-correlation function of the respective samples are constant as 
a function of scales, as expected from linear theory. All correlation 



functions appear to have a form that is well fit by a broken power- 
law with the dust-corrected distribution having a more pronounce 
convexity. It is noteworthy that the flux-limited analysis (Figure 
[5]l is more noisy than the volume-limited analysis (Figure |6), de- 
spite having a factor of two greater numbers. We believe that this 
is due to the inclusion of uncorrelated (in redshift space) galaxies, 
which dilute the clustering signal, especially for blue galaxies on 
small scales. From here onwards, we will restrict our analysis to 
the volume-limited sample. 

When we fit the projected correlation function of the full dis- 
tribution with the standard two parameter power-law for scales 
from 1 < rp < 10h~^ Mpc, red and green galaxies appears to 
have similar slope with 7 1.93, while blue galaxies have a 
substantially shallower slope, with 7 ~ 1.75. The best fit corre- 
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Figure 3. Contours and normalized counts of the two-dimensional correlation function C(7r, Tp) for red (left), green (middle) and blue (left) populations for the 
flux-limited sample. The panels on the top row are for the full distribution, bottom row are restricted to face-on galaxies (b/a > 0.6). The contours obtained 
after 2x2 Mpc boxcar smoothing. The levels are 0.0 (dotted lines), 0.25, 0.5, 1.0 (heavy Unes), 2.0, 4.0, 8.0 and 16.0. Red galaxies have the strongest 
finger of God effect (extension in the tt direction). The finger-of-God effect is seen in the green sample as well, while it is not present in the blue. Galaxies in 
the green valley appear to have clustering characteristics of the red sequence. 



lation lengths ro are 7.5h~^ Mpc (red), 5.3h~^ Mpc (green), and 
3.9h~^ Mpc (blue). The result of the covariance analysis using 999 
bootstrap samples is shown as confidence interval contours - 68% 
(inner) and 95% (outer) in F igure [7] (left p anel). These results are 
in excellent agreement with Heinis et al. I ( I2OO9 I). The correlation 
function from the dust-corrected distribution is substantially nois- 
ier, in part due to the smaller number of galaxies. Applying the 
same power-law fit, we obtain a larger uncertainty in 7, but compa- 
rable uncertainty in ro. This is shown on the right panel of Figure 
|7]and tabulated in Table 1. 

On the left panel, the projected correlation function of green 
galaxies is evidently intermediate between the red and the blue on 
a range of scales (0.2 < rp < 6h^^ Mpc), and runs faithfully 
parallel with the red correlation function, but with a lower ampli- 
tude. The picture for the dust-corrected distribution is slightly dif- 
ferent. While the green still runs in between the red and blue for 
Tp < lh~^ Mpc and mostly parallel to the red with a lower ampli- 
tude for broad range of scales, for larger scales (rp > 2h~^ Mpc), 
it coincides with the blue. This is a consequence of the more promi- 
nent two-halo excess on large scales seen in the blue correlation 
function. 

One way to compare the clustering bet ween subpopulation s 
of galaxies is through their relative bias (e.g. iNorberg et aT]|2002h . 



defined as 

where is the correlation function of interest, ^, is the fidu- 
cial correlation function that all correlation function is compare 
with. The approximate equa lity holds when the correlation func- 
tion takes a power-law form. lZehavi et alj ( |2005|) found that a typ- 
ical L* galaxies in the local universe with —20 < Mr < —21 
from the SDSS main spectroscopic sample have a fiducial power- 
law correlation function with ^, (r) = {r/5h^^ Mpc)^^'* on scale 
1 < rp < 10h~^ Mpc. By fitting the correlation function us- 
ing a constant 7 = 1.8 on scales 1 < rp < 10h~^ Mpc, we 
obtained the bias relative at rp = 5h~^ Mpc to a fiducial L* 
galaxies for each of our subpopulations: fered/fe* = 1-53 ± 0.08, 
bgreen/b, = 1.08 ± 0.06 and &biue/b, = 0.81 ± 0.06 for the full- 
sample, and fered/&. = 1.47 ± 0.10, bgreen/b, = 1.13 ± 0.09 and 
&biue/fe* ~ 0.92±0.09 for the dust-corrected samples (Table 1). In 
FigurelH we plot bi/K{r = 5h~^ Mpc) = {ro/5h'^ Mpc)^ as 
a function of the co-movi ng number density a nd compare our val- 
ues with those obtained bv lZehavietai] ( l2005h . Co-moving number 
density of our sample is estimated using the 1/Vmax method an d 
further corrected to match those obtained by IWvder et~aL I ( l2007b . 
While all the samples by construction have almost identical co- 
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Figure 4. Contours and normalized counts of the two-dimensional correlation function S;{tt, rp) for red (left), green (middle) and blue (left) populations for 
the volume-limited sample. The panels on the top row are for the full distribution, bottom panels are restricted to face-on galaxies (b/a > 0.6). The contours 
obtained after 2x2 Mpc boxcar smoothing. The levels are 0.0 (dotted lines), 0.25, 0.5, 1.0 (heavy lines), 2.0, 4.0, 8.0 and 16.0. Red galaxies have the 
strongest finger of God effect (extension in the n direction). Similar to Figure|3] but with lower signal-to-noise, the finger-of-God effect is seen in the green 
sample as well, while it is not present in the blue. 




Figure 5. The projected correlation function, uip(rp), of red sequence, green valley, and blue sequence galaxies for the flux-limited sample. The left (right) 
panel is from the full (dust-corrected) distribution. From red to blue, there is an increase in convexity, with the progressive shallowing of the slope on large 
scales, indicating an excess of two-halo signals at scale (rp > lh~^ Mpc). 
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Figure 6. The projected congelation function, Wp{rp), of red sequence, green valley, and blue sequence galaxies for the volume-limited sample. Similar 
to Figure \5\ from red to blue, there is an increase in convexity, with the progressive shallowing if the slope on large scale, indicating an excess of two- 
halo signals, at scale (rp > lh~^ Mpc). For the full distribution (left panel), the green Wp{rp) is parallel to the red with a lower amplitude for scales 
0.2 < Tp < 10h~^ Mpc. For the dust-corrected distribution (right panel), the blue Wp{rp) approaches the green on large scales. 
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Figure 7. Confidence contours of the power-law fits to Wp{rp) on scales Ih^^ Mpc < Vp < 10h~^ Mpc for the full (left panel) and dust-corrected (right) 
distributions from the volume-limited sample. The three colors are for red, green and blue subsamples, with the inner (outer) contours encircling the 68 (95) 
% confidence region in the ro-7 space. 



moving density, there is a range of relative bias, with the red hav- 
ing the highest bias, and blue the lowest. The spread among the bias 
are much smaller among the dust-corrected samples. Red galaxies 
have clustering strength above the nominal strength (inferred from 
a typical SDSS galaxy) for a given number density. One plausible 
scenario suggest that red galaxies have a larger than average satel- 
lite fraction. On the other hand, blue and green galaxies have lower 
observed number density compare to their clustering. If we assume 
all blue and green galaxies are central galaxies in the halo they re- 



side in, the lower observed number densities (~ O.OOl/i"^ Mpc~^) 
suggest that only a small fraction of these halos, as described by 
their lower bias, host a blue or green galaxy, as their expected num- 
ber density (inferred from SDSS to be ~ 0.01^'^ Mpc^'^) is much 
higher. These fractions are increased if the average phases of green 
and blue is shorter then the life time of the halos JHaiman & Huil 
l200lUMartini & Weinberj200lh . In the case of the green galaxies, 
the transitional nature of these galaxies may be closely related to 
the AGN they host, or to minor mergers and starbursts, with trig- 
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Table 1. Power-law fit to the projected auto-correlation function of volume-limited samples (Ih ^ Mpc < Vp < lOh ^ Mpc) 



Full-Sample N z median Mr ro 7 Relative Bias x^/dof 



Red 1971 (0.03,0.12) -20.9 7.5 ± 0.27 1.94 ±0.07 1.53 ± 0.08 1.36 

Green 1971 (0.03,0.12) -20.9 5.3 ±0.19 1.93 ± 0.08 1.08 ± 0.06 0.40 

Blue 1971 (0.03,0.12) -20.9 3.9 ±0.25 1.73 ±0.10 0.81 ± 0.06 1.54 



Face-On 

Red 1226 (0.03,0.12) -20.9 7.6 ± 0.38 1.81 ±0.11 1.47 ±0.10 2.47 

Green 1226 (0.03,0.12) -20.9 5.5 ±0.24 1.92 ±0.15 1.13 ±0.09 0.48 

Blue 1226 (0.03,0.12) -20.9 4.6 ± 0.31 1.74 ±0.17 0.92 ± 0.09 1.17 



4.5 Cross-Correlation Functions 
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Figure 8. Bias relative to the typical L* galaxy with a ^^"^(r) = 
(r/ro)~^ * at rp = 5h~^ Mpc versus co-moving number density. We 
also plotted points from SDSS (Zehavi et al.2005). The soHd circles are 
from their luminosity range analysis, i.e. galaxy with absolute magnitude 
between Mi and A/2, while the open circles are from their luminosity 
threshold analysis, i.e. Mr < Mu^n- While far from unique, the higher 
red bias suggests that on average red galaxies are more likely to be satellite 
galaxies then a typical galaxy with those number densities. If we assume 
the lower blue and green galaxies are all central galaxies in the halo they 
reside in, the lower observed number densities suggest that only a fraction 
of these halos, as described by their bias, host a blue or green galaxy, as 
their expected number density is much higher. 



gering cycles corresponding to the duty cycles for these respective 
phenomena. We defer a more detail analysis of halo occupation and 
the life cycle of transition galaxies using star-formation rate tracers 
of different lags to a future paper (Heinis et al. in prep.). Using 
a theo retical bias function of dark matter halos jSeliak & WarrenI 
I2OO4I) and normalizing to A/* = lO^^ft"^ MgB red galaxies clus- 
ter similar to halos with mass 



10^2-^/i~^ Mq, and blue galaxies 10 



IQ^'^ 'h ^ Mq, green galaxies 



The halo mass that a typical L* galaxy resides in. 



Figure |9] shows the two-dimensional C(rp,7r) cross-correlation 
function (CCF) for the three cross-pairs of galaxies in our volume- 
limited sample, with the panels on the top for the full distribu- 
tion, and the bottom panels for the dust-corrected distribution. The 
large scale infall effect is seen in all three panels suggesting that 
the galaxies, on average, trace a similar matter distribution as ex- 
pected from linear theory. On small scales, the finger-of-God ef- 
fect is strongest for the red-green CCF, intermediate for red-blue 
and green-blue. In Figure [TO] we plot the relative amplitude of ^ 
at the projected distance Vp — lh~^ Mpc as a function of line- 
of-sight distance tt. The circles are for green-red cross correlation, 
squares for red-blue, and triangles for green-blue. Green-red have 
largest amplitude for a wide range of tt, while red-blue and green- 
blue have smaller amplitudes. The stronger finger-of-god compo- 
nent suggest that the dynamics of red and green are more strongly 
coupled compared with red-blue or green-blue pairs. Note that red- 
blue CCF have a very different ^(tt, Vp) compared with the green 
auto-correlation function (ACF), the former has a weaker finger-of- 
God and stronger infall compression. 

As discussed in section [3. 1.21 one way to understand the re- 
lationship between two populations is to compare the projected 
cross-correlation function between the two with their geometric 
mean. If the two populations are mixed evenly, the cross-correlation 
functions should trace the geometric mean. If they are spatial seg- 
regated (partially) beyond what was expected from their respec- 
tive auto-correlation function, the cross-correlation function should 
be systematically below the geometric mean. The projected cross- 
correlation functions for the red and blue galaxies Wp^{rp) are 
plotted as open circles in Figure [TT] Also plotted for comparison 
are the red and blue auto-correlation functions (solid lines), and 
the red-blue geometric mean iUp''(rp) (dashed lines). On scales 
larger than ~ 'ih^^ Mpc, the Wp'^{rp) approaches the geometric 
mean. For scales 0.2 < rp < 2h~^ Mpc, Wp'^ (r-p) is systemati- 
cally below the geometric mean for the full distribution (left panel). 
For the dust-corrected distribution (right panel), ''(rp) starts to 
inch below the geometric mean only from Vp < 0.7/i~^ Mpc on- 
wards. Our results are consistent with the partia l morphology seg- 
regation within galaxy clusters l lDressleilll980l) . At the one-halo 
regime (rp < 1 — 3/i~^ Mpc), the relevant scales for galaxy clus- 
ters, red galaxies tend to occupy the cores of the clusters, while 
blue galaxies tend to lie towards the periphery. The lower level of 
spatial mixing on these scales suppresses amplitude of the cross- 
correlation function. Note that the green auto-correlation function 
(green solid line) does trace the red-blue geometric mean on large 
scales, suggesting that blue and red galaxies do mix evenly on these 
scales. 
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Figure 9. Contours and normalized counts of the two-dimensional cross-correlation function §(71, rp) for red-blue (left), green-red (middle) and green-blue 
(left) for the volume-limited sample. The panels on the top row are for the full distribution, bottom panels are restricted to face-on galaxies {b/a > 0.6). The 
contours obtained after 2x2 h^^ Mpc boxcar smoothing. The levels are 0.0 (dotted lines), 0.25, 0.5, 1.0 (heavy lines), 2.0, 4.0, 8.0 and 16.0. The large scale 
infall effect is seen in all three panels. On small scales, the finger-of-God effect is strongest for the red-green, intermediate for red-blue and weakest for green- 
blue, as expected. Note that red-blue CCF ^(Tr, rp) has a weaker finger-of-God and stronger infall compression compare to the green valley auto-correlation 
function (Figure|4]middle panels). 



Figure [T2I shows the projected cross-correlation function be- 
tween the green and red w|'' (r^), and green and blue Wp'^{rp). 
The solid green line is the auto-correlation function of the green 
sample {Wp^[rp)), while the dashed lines are the geometric means: 

'^p'^iXv) '^p'^i'^'p)- Fo'' the full distribution (left panel), 
'Wp''^{rp) is systematically below Wp'"^ for a range of scales within 
errors. ui|''(rp) are consistent, or slightly above Wp'''{rp) for rp > 
lh~^ Mpc, and systematically below for 0.3 < rp < Ih"^ Mpc. 
These inferences are shown more clearly in Figure [T3] (left for full, 
right for dust-corrected), where we plot the normalized cross clus- 
tering strength - the ratio Wp'^ /wp^ , where x is either red, or blue 
- relative to the auto-correlation of the green galaxies. The addi- 
tional solid lines are the normalized ACF of the red and blue. This 
suggest that red and green are consistent with being drawn from the 
same statistical sample on average for a large range of scales with 
a slight anti-bias on small scales. For the green and blue, there is 
a stronger anti-bias for r < lh~^ Mpc as they avoid each other 
on these scales. Our results suggest that on scales typical of dark 
matter halos, galaxies drawn from the green and blue populations 
are less associated than would be predicted by their respective auto- 
correlation function. 



5 DISCUSSION 

5.1 Comparison with DEEP2 2 ~ 1 

The GALEX NUV selected sample used in our analysis is di- 
rectly comparable t o the high redshift (z ~ 1) galaxy sample of 
the DEEP2 survey dCoil et alJlioOS) since their optical selection 
mimics the rest-frame NUV. In that paper, clustering analysis was 
done on green valley galaxies for the first time. In their analysis, 
the green valley S,{rp,n) appears to display kinematic structure 
intermediate between red and blue galaxies, with an intermediate 
finger-of-God effect and an intermediate overall clustering ampli- 
tude. When redshift space distortion is removed, the projected cor- 
relation function Wp{rp) shows a scale dependence convergence. 
At rp > \h^^ Mpc the functions converge to those of red popula- 
tions, while for rp < Ih^^ Mpc it tends toward the clustering of 
the blue population. 

In contrast, we find that our green valley £,{rp, n) has a strong 
finger-of-God effect consistent with that measured for red sequence 
galaxies, differing only in their amplitude. This can be seen most 
clearly in the Wp{rp) analysis for the volume-limited sample (Fig- 
ure |6j, where the green function has form similar to that of the red 
function for scales 0.2 < rp < 15h~^ Mpc, but is displaced to 
lower amplitude. It is noteworthy that the slope of the blue cor- 
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Figure 11. The projected cross-con'elation function between red and blue sequence galaxies Wp'^(rp) are shown as open circles. The solid lines are the 

auto-coiTelation function of the red and blue galaxies. The dashed line is the red-blue geometric mean uip ''(r^). On scales larger than ~ 3h~^ Mpc, the 
cross-correlation approaches the geometric mean. For scales 0.2h~^ Mpc < rp < 2h~^ Mpc. the cross-correlation function is systematically below the 
geometric mean. This result suggests that on small scales, blue an d red galaxies a re spatially segregated beyond what was expected from their auto-correlation 
function. This is consistent with the morphology density relation )DressleJl980h . 




0.1 1.0 10.0 0.1 1.0 10. 
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Figure 12. The projected cross-correlation function of green sample with blue uip'''(r'p) (blue open circles) and red Wp'^ ){fp) (red open circles) subsamples. 
The solid line shows the auto-correlation function of the green sample. The dashed lines are the geometric mean of the green-red (in red) and green-blue (in 
blue) auto-correlation. On the left (right) we show the analysis from the full (dust-corrected) distribution from the volume-limited sample. (See Figure [TSl for 
more details.) 



relation function begins shallowing at r ~ lh~^ Mpc, display- 
ing the kind of one-halo and two-halo segregation expected from 
a correlation fun ction dominated by central galaxies. In contrast to 
et al.l jlOOSh . the green auto-correlation function converges to 
that of the blue at Vp > lh~^ Mp c. We emphasize that the green 
galaxy sample of C oil et aL||200g) is defined differently from ours. 
We divide the CMD into three disjoint parts to separate our galax- 



ies into the subpopulations while in et al.L the green overlaps 
both the red and the blue. 



5.2 Green Valley 

Many recent studies reveal that blue sequence mass has remained 
roughly constant since 2 ~ 1 telantorj|2006l ; iFaber et alj|2007h 
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Figure 13. The ratio of projected congelation function {Wp{rp)) with the auto- congelation function {Wp^ {vp)) of the green subsample. The two ratios from 
the cross-conelation functions, Wp'^ [rp] and Wp' (vp), are indicated by colored open circles. Dashed lines are for geometric means, as in figure [T2I The 
additional solid lines are ratios from the auto-correlation function of red and blue respectively. The fact that both Wp''^ (rp) and Wp'^{rp) are below their 
respective geometric mean on small scales suggest that both red and blue galaxies are found in different environments compared to the green on these scales, 
albeit with large uncertainties. 
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Figure 10. This figure shows the relative amplitude of ^ at projected dis- 
tance Tp = lh~^ Mpc as a function of the line-of-sight distance tt. The 
circles are for green-red cross correlation, squares for red-blue, and trian- 
gles for green-blue. Green-red have largest amplitude for a wide range of tt, 
while red-blue and green-blue have smaller amplitudes. These results sug- 
gest that the finger-of-god component of green-red cross correlation func- 
tion is stronger then those from red-blue or green-blue. 



because the average ongoing star formation over < z < 1 is 
balanced by mass flux off the blue sequence, pre sumably towards 
the build-up of th e red sequence since z ~ 1 teelletaLlllOoi : 
iMartin et alj2007h . Hence, green valley galaxies occupy a position 



where one expect to find many transitional galaxies. iMartin et al.l 
also note that the AGN fraction peaks at the green valley. For our 
green valley definition using the oblique color cuts from Figure [T] 
the AGN fraction is ~ 509Cl 

In their study of AGN using SDSS, IConstantin & VogelevI 

( l2006l) found that the redshift space two-point correlation function 
of Seyferts is less clustered than that of LINERS (low-ioniz ation 
nuclear emissio n-line regions). However, iMiller et al.l ( l2003h and 
iLi et alj ( l2006bl) found that AGN as a whole cluster similarly as 
typical L* galaxies , if on e takes into account luminosity bias. 
l\Vang & KauffmannlJiooi) argues that almost all galaxies in the lo- 
cal universe with stellar mass > 10h~^ Mq have active nuclei, of- 
ten LINERS with lines too weak to be detected spectroscopically in 
SDSS. Our green galaxies, with luminosities peaked at Mr ~ —21, 
have stellar masses well above 10h~^ could be domi- 

nated by LINERS (either detected or undetected). This would in 
part explain the clustering effect we see in Figures[3]and|4](middle 
panels) where those bulge-dominated LINERS display kinematics 
similar to red sequence (primarily non-active) gala xies. The rea- 
son that AGN from r-band selected survey (e.g. Mill er et al.ll2003h 
cluster on average much like typical L* galaxies may merely be 

coincidental. 

The lack of such behavior in the z ~ 1 green valley jCoil et al] 

I2OO8I) may be attributable to evolution in the AGN population. 
There may be fewer LINERS, or the red sequence galaxies may 
have been experiencing more gas infall, feeding their AGN. We 
note that with a relative bias ~ 1.1, green galaxies cluster similarly 
to a typical galaxy with L* — 0.5 ~ —21, the median luminosity 
of our sample. 

As was discussed in section HTl a substantial fraction of the 



^ This is a lower limit since the fraction of low luminosity AGN and com- 
posite objects is unknown. 



UV- Optical Clustering 15 



galaxies in the green valley are dusty star forming galaxies. Be- 
cause dust content can modify the color-magnitude diagram we 
use to separate the galaxies into red, green and blue populations, 
this might potentially change the behavior of the correlation func- 
tion. We argue here that to the extent that dust modifies the number 
counts of green galaxies, it is to pr omote the migra tion from the 
blue sequence to the green valley IChoietal.ll2007h . Our results 
from Figure |4] (for the volume-limited sample) suggest that their 
influence in modest at best, and merely acts as additional poisson 
noise to the two-dimensional ^(tt, rp) without altering the kinemat- 
ics in the sample. The projected correlation functions of Figure |6] 
show similar characteristics. 

5.3 Blue (UV) vs. Blue (Optical) 

In their studies using g — r color. lZehavi et alj (^005') found that the 
correlation function of blue galaxies exhibits a lower amplitude and 
shallower correlation functions. By fitting the Wpij-p) with halo oc- 
cupation distribution (HOD) models, they found that the majority 
(~ 70 — 90%) of blue galaxies are central galaxies in dark mat- 
ter halos, usually halos with mass < 10^^ Mq, in contrast to only 
the most luminous red galaxies being central objects of massive 
halos, while the majority of (less luminous) red galaxies are satel- 
lites. This fits well with the notion that blue galaxies are field galax- 
ies - the central objects of low mass halos; red galaxies, with the 
exception of the central galaxy in clusters and groups, are mostly 
satellites of massive halos. 

One can infer a similar conclusion from the auto-correlation 
Wp{rp) plot for the blue galaxies (Figure|6j. If we decompose the 
correlation function into two parts, one due to the one halo term and 
the other due to the two halo term, blue galaxies show strong two- 
halo excess on scales > lh~^ Mpc, implying that on those such 
scales, the majority of galaxy pair s are from different halos. This 
was seen in the optical analysis of IZehavi et alj jlOOSt) but promi- 
nently in our blue sequence sample. For the dust-corrected volume- 
limited sample (right panel of Figure|6j, the amplitude of the auto- 
correlation function of the blue sample beyond ~ Sh~^ Mpc actu- 
ally rises to match th e clustering strength of the green sample. 

iLi et alj l l2006ah found that the dependence of Wp{rp) on op- 
tical g — r color extends beyond 5/i~^ Mpc, suggesting that the 
conventional wisdom that clustering should converge at large scales 
may not occur until at a scales larger than 5h~^ Mpc. Here, we ar- 
gue that this is an effect due to the mixture of population between 
blue galaxies, and green and red; that the optical g — r color does 
not have sufficient power to separate the green from the blue. Blue 
galaxies, by themselves, have a very pronounced two halo excess 
and are dominated entirely by central galaxies, flattening the cor- 
relation function substantially at large scales, compared to the red 
and green population. To the extent that one can eliminate or cor- 
rect for the internal reddening due to dust, NUV — r color is very 
efficient in isolating a sequence of purely star-forming galaxies. 



6 SUMMARY AND CONCLUSIONS 

We have constructed a GALEXand SDSS matched catalog, where 
we have used the GR3 catalog from GALEXand the SDSS DR5 
main spectroscopic galaxy sample. We construct the galaxy distri- 
bution of NUV — r vs Air color-magnitude diagram, and divide 
the distribution into populations of red sequence, green valley and 
blue sequence. Since our main goal is to study the color dependence 
of clustering, we took substantial care in matching the luminosity 



distribution of each population. For each population, we measure 
the two-dimensional correlation function ^(Tr, rp), and the one- 
dimensional projected correlation function Wp. We also perform 
cross-correlation analyses between each of the sub-populations. 

Our principle finding is that the red sequence and green val- 
ley appear to show similar clustering properties, as expressed in the 
finger-of-God effect in the auto-correlation function. The projected 
correlation function is consistent with red and green galaxies re- 
siding as satellites of massive halos, while the blue sequence shows 
what appears to be a clear two-halo signature, hence primarily serv- 
ing as central galaxies of less massive halos. The cross-correlation 
function also shows that green and blue galaxies, on small scales, 
are not a mere statistical mix, but are spatially segregated from each 
other. 

The findings would appear to place the green valley population 
with the red sequence. The green valley would largely consist of 
massive galaxies that reside in massive halos, and whic h cluster like 
the re d seq uence. We note tha t iMartin et alj J2007l) . IWvder et al.l 
I2OO7I), and ISalim et al] l l2007h show that a large fraction of type 
II AGN are found in the green valley. ISalim et aU show that in the 
plot of specific star formation rate vs stellar mass, the AGN tend to 
be found in massive (> 3 x 10^" Af©) galaxies. The AGN occupy 
a region in these plots that strongly resembles that of the reddest 
class of galaxies, the "no-Ha" red sequence galaxies. Significantly, 
the AGN are clearly offset from the locus of the blue sequence, in 
the plot of specific star-formation rate (SFR) versus mass. The sig- 
nificance is that while a minority of AGN are found with properties 
that coincide with those of the more massive blue sequence galax- 
ies, green valley galaxies - the subsample with the largest AGN 
fraction — exhibit properties similar to those of the red sequence, 
but showing mildly elevated star formation. 

In this study, we have shown that the green valley population 
clusters in ways that are characteristic of, but also less strongly 
than, the red sequence. One may suggest that these studies paint 
a picture in which both the properties of the green valley and the 
"demographics" are different from those of the blue sequence, at 
the present epoch. These findings do not necessarily contradict the 
studies that find an increase in the total mass of red sequence galax- 
ies since z ^ 1. They do suggest, however, that if blue sequence 
galaxies evolve by some process to the green valley, and ultimately 
to the red sequence, such evolution must also accompanied by a 
transition from the field environment to a group/cluster environ- 
ment. Such a change could conceivably occur if the blue popula- 
tion resides along filaments that infall into clusters, over time. Our 
cross-correlation results show that green galaxies avoid both red 
and blue galaxies on small scales is consistent with the change 
in environmen t hypothesis. We no te that models like ram pres- 
sure stripping jGunn & Gottlll972h, starvation of (cold ) gas, and 
the virial shock heating model of ( Dekel & Birnboim 2006) nat- 
urally incorporate environmental factors in their mechanism for 
color transformation in galaxies. 

It is also possible that the downsizing i lCowie etalJl 19961) ef- 
fect is so strong th at most star fo rming galaxies are evolving rapidly 
with redshift (e.g. lTinslevlll968h . One must recall that star forming 
activity at 2; ~ 1 resides in considerably more massive galaxies, 
and that a color-based population separation, as we have done, will 
refer to much higher masses; the rest-frame colors may be similar, 
but the fundamental nature of the galaxies, not. 

One may speculate that the the green valley is occupied by 
nominally red galaxies that experience the infall of a gas rich sys- 
tem that either induces star formation and/or fuels the AGN, ren- 
dering it visible via its emission lines. However, the feedback of an 
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AGN might inhibit star formation and move a blue sequence galaxy 
to the green valley. Any number of environmental effects (e.g. ha- 
rassment, starvation) might speed the consumption of gas in a disk, 
again moving a galaxy to the green valley. A small sample of opti- 
cally quiescent members of the green valley that nominally have a 
UV excess show clear star f ormation signatures (spirals) when im- 
aged in the UV using HST jSalim & Richll2009h . There is still the 
issue of the origin of the low mass red sequence, and the evolution 
of blue sequence, into the green valley and ultimately the red se- 
quence, might have an important role in the growth of the lower 
mass portion of the red se quence. This was p artial ly addressed 
by sem i-analytical work of iBenson et al.l ( |2003|) and iBower et al.l 
( I2OO6I) . 

In contrast to the construction of color-magnitude diagrams 
for stellar populations, the environment, for galaxies, is a critical 
physical variable in their evolution. This is true both the in the 
sense of their dark matter environment as well as the presence of 
detectable companion stellar systems. In considering the major pro- 
cesses driving galaxy evolution, it would appear that evolution of 
both of these observables must be considered, as a function of look- 
back time. 



ACKNOWLEDGMENTS 

Y.S.L. would like to thank C. Hirata, S. Salim, C. Park, J. Kor- 
mendy and Z. Zheng for helpful discussions. This work has made 
extensive use of IDLUTIL^ and Goddard IDL libraries. RMR ac- 
knowledges support from grant GO-11182 from the Space Tele- 
scope Science Institute. 

GALEX(Galaxy Evolution Explorer) is a NASA Small Ex- 
plorer, launched in April 2003. We gratefully acknowledge NASA's 
support for construction, operation, and science analysis for the 
GALEXmission, developed in cooperation with the Centre Na- 
tional dEtudes Spatiales of France and the Korean Ministry of Sci- 
ence and Technology. 

Facilities: GALEX, SDSS 



REFERENCES 

Basu-Zych, A.R., et al., 2008, submitted to ApJ 

Bell, E.R, et al., 2004, ApJ, 608, 752 

Benson, A.J., et al., 2003, ApJ, 599, 38 

Bertin, E., & Arnouts, S., 1996, A&AS, 117, 393 

Blanton, M.R., 2000, ApJ, 544, 63 

Blanton, M.R., et al., 2003, ApJ, 592, 819 

Blanton, M.R., et al., 2005, AJ, 129, 2562 

Blanton, M.R., et al., 2005, ApJ, 629, 143 

Blanton, M.R., 2006, ApJ, 648, 268 

Blanton, M.R., & Roweis, S., 2007, AJ, 133, 734 

Bower, R.G., et al., 1992, MNRAS, 254, 601 

Bower, R.G., et al., 2006, MNRAS, 370, 645 

Brown, M.J.I., et al., 2008, ApJ, 682, 937 

Budavari, T, et al., 2003, ApJ, 595, 59 

Chen, J., 2009, Ac&A, 494, 867 

Choi, Y.-Y., et al., 2007, 658, 884 

Coil, A.L., et al, 2008, ApJ, 672, 153 

Constantin, A., & Vogeley, M.S., 2006, ApJ, 650, 727 

Cowie, L.L, et al., 1996, AJ, 112, 839 

* http://spectro.princeton.edu/idlutils 



Croton, D.J., et al., 2006, MNRAS, 365, 11 

Davis, M., & Geller, M.J., 1976, ApJ, 208, 13 

Davis, M., & Peebles, PJ.E., 1983, ApJ, 267, 465 

Davison, A.C., & Hinkley, D.V., 1997, Bootstrap Methods and 

Their Applications (Cambridge: Cambridge Univ. Press) 
Dekel, A., & Birnboim, Y, 2006, MNRAS, 368, 2 
Dressier, A., 1980, ApJ, 236, 351 
Efron, B., 1981, Biometrika, 68, 589 
Efstathiou, G., et al., 1991, ApJ, 380L, 47 
Faber, S.M., et al., 2007, ApJ, 665, 265 
Gunn, J.E., & Gott, J.R., 1972, ApJ, 176, 1 
Heinis, S., et al., 2007, ApJS, 173, 503 
Heinis, S., et al., 2009, submitted to ApJ 
Haiman, Z., & Hui, L., 2001, ApJ, 547, 27 
Hamilton, A.J.S., 1992, ApJ, 385, L5 
Hamilton, A.J.S., & Tegmark, M., 2002, MNRAS, 330, 506 
Hopkins, PR, et al., 2007, ApJS, 163, 1 
Hopkins, PR, et al., 2007, ApJ, 659, 976 

Hubble E.P et al., 1936, The Realm of the Nebulae (Oxford: Ox- 
ford Univ. Press), 79 
Johnson, B.D., et al., 2006, ApJ, 619, L109 
Kauffmann, G. et al., 2004, MNRAS, 341, 33 
Kaiser, N., 1987, MNRAS, 227, 1 
Kerscher, M., et al., 2000, ApJ, 535, L5 
Kron, R. G. 1980, ApJS, 43, 305 
Landy, S.D., & Szalay, A., 1993, ApJ, 412, 64 
Li, C, et al., 2006a, MNRAS, 368, 21 
Li, C, et al., 2006b, MNRAS, 373, 457 
Loh, J.M., 2008, ApJ, 681, 726 
Madgwick, D.S., et al., 2003, MNRAS, 344, 847 
Martin, D.C., et al., 2005, ApJ, 619, LI 
Martin, D.C., et al., 2007, ApJS, 173, 342 
Martini, P, & Weinberg, D.H., 2001, ApJ, 547, 12 
Miller, C.J., et al., 2003, ApJ, 597, 142 
Milhard, B., et al., 2007, ApJS, 173, 
Morrissey, P, et al., 2005, 619, L7 
Morrissey, P, et al., 2007, ApJS, 173, 682 
Norberg, P, et al., 2002, MNRAS, 332, 837 
Norberg, P, et al., 2009, arXiv:0810.1885vl 
Ostriker, J.P, & Steinhardt, PJ., 1995, Nature, 377, 600 
Padmanabhan, N., et al., 2008, arXiv :O80272TO5l 
Park, C, et al., 2007, ApJ, 658, 898 

Peebles, P.J.E., 1980, The Large-scale Structure of the Universe, 

Princeton Univ. Press 
Perlmutter, S., et al., 1999, ApJ, 517, 565 
Rich, R.M., et al., 2005, ApJ, 619L, 107 
Salim, S., & Rich, R.M., 2009, BAAS 213, 435.03 
Riess, A.G., et al., 1998, AJ, 116, 1009 
Sahm, S., et al., 2007, ApJS, 173, 267 
Sargent, W.L.W., & Turner, E.L., 1977, ApJ, 212, L3 
Seljak, U., & Warren, M.S., 2004, MNRAS, 355, 129 
Schinomovich, D., et al., 2007, ApJS, 173, 315 
Silk, J., & Rees, M.J., 1998, A&A, 33 IL, 1 
Swanson, M.E.C., et al., (2008), MNRAS, 387, 1391 
Spergel, D.N., et al., 2007, ApJS, 170, 377 
Tinker, J., 2007, MNRAS, 374, 477 
Tinsley, B.M., 1968, ApJ, 151, 547 
Toomre, A., & Toomre, J., 1972, ApJ, 178, 623 
Totsuji, H., & Kihara, T, 1969, PASJ, 21, 221 
Wang, L., & Kauffmann, G., 2008, MNRAS, 391 785 
Wang, Y, et al., 2007, ApJ, 664, 609 
Wyder, T. K., et al.2007, ApJS, 173, 293 



York D. et al., 2000, AJ, 120, 1579 
Zehavi, I., et al., 2005, ApJ, 630, 1 

Zwicky, F., et al., 1968, Catalog of Galaxies and of Clusters of 
Galaxies, Vols. 1-6 (Pasadena, Caltech) 



UV- Optical Clustering 



